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ABSTRACT 



Histogram-reweighting Monte Carlo simulations were used to obtain polymer / sol- 
vent phase diagrams for lattice homopolymers of chain lengths up to f=1000 monomers. 
The simulation technique was based on performing a series of grand canonical Monte 
Carlo calculations for a small number of state points and combining the results to obtain 
the phase behavior of a system over a range of temperatures and densities. Critical pa- 
rameters were determined from mixed-field finite-size scaling concepts by matching the 
order parameter distribution near the critical point to the distribution for the three- 
dimensional Ising universality class. Calculations for the simple cubic lattice 
(coordination number z=6) and for a high coordination niraiber version of the same lattice 
(z=26) were performed for chain lengths significantly longer than in previous simulation 
studies. The critical temperature was found to scale with chain length following the 
Hory-Huggins functional form. For the z=6 lattice, the extrapolated infinite chain length 
critical temperature is 3.70±0.01, in excellent agreement with previous calculations of the 
temperature at which the osmotic second virial coefficient is zero and the mean end-lo- 
end distance proportional to the number of bonds. This confirms that the three alternative 
definitions of the © temperature are equivalent in the limit of long chains. The critical 
volume fraction scales with chain length with an exponent equal to 0.38+0.01, in agree- 
ment with experimental data but in disagreement with polymer solution theories. The 
width of the coexistence curve prefactor was tentatively found to scale with chain length 
with an exponent of 0.20±0.03 for z = 6 and 0.22±0.03 for z = 26. These values are near 
the lower range of values obtained from experimental data. 
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INTRODUCTION 



Phase equilibria in polymer solutions are important in manufacturing, processing and 
applications of macromolecules. Significant progress has been made in recent years in 
the development of detailed atomistic models that can be used to predict properties of 
polymeric systems. Phase coexistence properties, however, are generally difficult to ob- 
tain directly for atomistically detailed models because the free energy of a system cannot 
be easily determined from simulations. Simple lattice models have often been used to 
obtain results that can be directly compared with statistical mechanical theories, but 
even for simple models, there is only a limited number of previous studies of coexistence 
properties for polymer/solvent systems. 

Results for phase coexistence and critical properties of relatively short lattice ho- 

mopolymers have been obtained previously by Yan et al. ' on the z=6 simple cubic lattice 
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for chains of length r up to 200, by Mackie et al. for the z=26 lattice for chains up to 

r=128 and by Wilding et al.^ for the bond fluctuation model for chains up to r=60. Co- 

4 

existence curves for continuous-space models were obtained by Sheng et al. for a bead- 
spring model for chains up to r=100 and Escobedo and de Pablo' for square- well chains 
up to (=100. These previous studies have confirmed that the critical temperature. 7^., 
depends on chain length r in a manner consistent with the functional form suggested by 
Flory-Huggins theory, 

^ (1) 

TJr) r,(~) V7 2r 

where r^.('=) is the critical temperature for chains of infinite length. For long chain 
lengths, equation 1 reduces to T^.(°°)-T^.(r) « r with jr3=0.5, following the notation of 
reference*. 
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Another important scaling relationship is that for the chain length dependence of the 
critical volume fraction, (|)^. , 

^,(r)^r-'\ (2) 

In previous lattice-based studies,' ' the exponent X2 was found to be near or below the 
experimentally determined values' ** of 0.38 to 0.40 and significantly lower than the value 

X2 = 1/2 suggested by Flory-Huggins theory. For the previous study of the continuous- 

4 

space bead-spring model the longest chain lengths studies were probably not in the 
scaling regime and the apparent exponent was significantly lower. 

Finally, the chain-length dependence of the width of the coexistence curve near 
the critical point is expected to be described by*, 

(r, T) - is,2 (r, T) « r-'' (l-T/T, (r))^ , (3) 
where and (^2 are the volume fractions in the two coexisting phases and p=0.326 is the 
universal critical point exponent appropriate for three-dimensional fluid systems. The 
exponent jci has been determined experimentally*"' to be between 0.23 and 0.34 but no 
simulation estimates are available to the best of our knowledge. 

In the present study, the histogram reweighting grand canonical Monte Carlo 

9 10 
simulation technique combined with mixed-field finite scaling concepts has been em- 
ployed. This technique has been recently applied by Wilding et al^ to calculate poly- 
mer/solvent critical point parameters of the bond fluctuation model of chains of length up 
to r=60. Because of the higher flexibility of the bond fluctuation model, these chain 
lengthss are equivalent to significantly longer ones on the simple cubic lattice. 

The first part of the present manuscript deals with methodological issues related to 
the application of histogram-reweighting grand canonical Monte Carlo simulations to the 
prediction of phase diagrams and critical points of lattice homopolymers. The following 
section presents results for cubic lattices of coordination number z=6 and z=26. Results 
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for the z=6 lattice are only in modest agreement with previous calculations for chain 
lengths r=100 and 200. The infinite-chain length critical temperature is compared to in- 
dependent estimates of the temperature at which chain dimensions behave ideally and the 
chain-chain second virial coefficient vanishes. We also obtain estimates of the exponents 
for scaling with chain length of the critical temperature, critical volume fraction and co- 
existence curve width and compare the results to available experimental data. 

SIMULATION METHODS 

Histogram-Reweighting Monte Carlo 

9 

The method has been described previously ; here, we would like to summarize the tech- 
nique as applied to systems of interest to the present study. A grand canonical Monte 
Carlo simulation is performed in a simulation cell of size V under periodic boundary con- 
ditions, at an imposed value of the chemical potential n and a temperature T. Particles are 
created and annihilated using the standai'd acceptance criteria. ' ' The frequency of occur- 
rence, /fA'.ii), of N particles with total configurational energy E in the simulation cell is 

/( N, E) = Q.( N,V,E)e^^^-'^^ / S(vl,VJ) , (4) 

where Q.(N,V,E) is the microcanonical partition function (density of states), p is the in- 
verse temperature (|3 = kgT , where is Boltzmann's constant), and H is the grand 
canonical partition function. Given the distribution function fiN.E), collected in histo- 
gram form in the production period of a simulation, an estimate of the ratio of micro- 
canonical partition functions for the system under study for two different values of N in 
the range covered by the simulation can be obtained directly as 

Sii{Ni,V,Ei) ^ /(iVi,£i) ^3[n(Ar,-Af,)-(E,-£,)i 
Q.2{N2,V,E2) f{N2,E2) 

In addition, one expects that a simulation at a different value of the chemical potential, 
|i' , and temperature, T' , would result in a new distribution function, / '{N, E) , with 
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f(N,E) 

The rescaling suggested by equation 6 can only be performed over a limited range of 
chemical potentials and temperatures since the original simulation provides statistically 
significant results only over a finite range of particle numbers and energies. For extend- 
ing the range of particle numbers over which the partition function ratio can be deter- 
mined from equation 5, several nms wiU need to be performed at different values of the 
chemical potential that result in overlapping distribution functions fiN,E). From equation 
4, the microcanonical partition function over the range of densities covered in each indi- 
vidual run, with index n, can be obtained from 

In a„{N,V, E) = ln[/(iV, £)] - puiV + p£ C„ (7) 

where C„ is a run-specific constant equal to the logarithm of the grand partition function 
for the chemical potential and temperature of runn, \n's.(\i„,V ,T„) . 

To obtain an estimate of the microcanonical partition function valid over a broad range of 

particle numbers and energies, results from different simulations need to be combined by 

assigning values of C„ for each run in a self-consistent fashion. For combining results 

12 

from multiple runs, the technique of Ferrenberg and S wendsen is used. The probability 
/'(W,£;H,P) of a certain number of particles and a certain energy resulting by combining 
nms n=l through R, assuming that they all have the same statistical efficiency is'^ 

P(iV,£;it,P)^ ,F-^"^^;^^--'''-P^"P^' (8) 
exp[- KE + K^„N - C„] 

where is the total number of observations for run m. The constants C„ are obtained 
from an iterative relationship: 

exp[c„] = II/'(iV,£;n„,p„) (9) 

E N 
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The Ferrenberg-Swedsen method ensures that there is minimum deviation between ob- 
served and predicted histograms from the combined runs. 

Mixed-Field Finite-Size Scaling 

3 

In order to obtain critical parameter estimates, mixed-field finite size scaling methods 
were used. A series of grand canonical simulations were performed near the expected 

critical point. The resulting histograms were combined according to equations 8 and 9 to 
obtain self-consistent estimates of the distribution functions P(A',£;)x,p). According to 
finite-size scaling theory, one needs to define an ordering operator, 7f(, combining the 
number of particles N and energy E, 

■M = N-sE (10) 

where * is a non-universal "field mixing" parameter controlling the strength of coupling 
between energy and density fluctuations near the critical point. At the critical point, the 
normalized probability distribution at a given system size L = v"^, /"^(jc), assumes a 
universal shape, with x = a{L,r)y.('}K - M^) . The non-universal scale factor a(L,r) is 

chosen to result in unit variance for the distribution Pi(x). An example of the matching 

13 

of some of our data to the universal curve obtained from is shown in Figure 1. There is 
excellent agreement between our data and the universal curve even though there are rela- 

3 

tively few (<100) chains in the system, as also observed previously . 
Configurational-Bl\s Sampling 

For long polymers at moderate and higher densities, it is not practical to perform inser- 
tions and removals directly, as the probability of inserting a molecule in a random internal 
configuration without overlapping with existing particles decays exponentially with chain 

length. We have used the simple "athermal" version of configurational-bias sampling 
14 15 

methods ' to enable insertions and removals of the chain molecules. According to this 
scheme, the first monomer of a chain to be inserted is placed in a random position; if the 

position is occupied the attempted insertion fails. Subsequent monomers are placed on 
unoccupied positions on the lattice, provided such positions exist along the directions of 
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growth permitted by lattice connectivity. The "Rosenbluth weight" for each growth step 
is calculated as the ratio of the number of unoccupied sites divided by the total number of 
sites, z=6 or 26 depending on the lattice. The Rosenbluth weight of the chain when it has 
been fully grown, Wncw, is the product of the weights calculated during each growth step. 
The attempted insertion is accepted with probability that takes into account the energy 
change for the insertion and the Rosenbluth weight". The reverse occurs during removal: 
a chain selected at random for a trial removal is "unzipped" from one randomly selected 
end to the other, and the Rosenbluth weight of the existing configuration is calculated as 
segment by segment is eliminated. We have chosen to use this limited form of configu- 
rational-bias sampling because it is computationally quite efficient. 

Models and Simulation Details 

In the model employed, homopolymers of chain length r exist on a cubic lattice of coor- 
dination number z=6 or z=26. Monomeric solvent particles fill all space not occupied by 
polymer segments. For the z=6 case, non-bonded monomers interact when they are 
within one lattice spacing along the principal directions of the lattice: relative position 
vectors for interactions are (1,0,0) and the five additional vectors that result from sym- 
metry operations in the -x, y, -y, z and -z directions. For the z=26 lattice, the relative po- 
sition vectors for interactions are (1,0,0) (1,1,0), (1,1,1) and additional vectors that result 
from symmetry operations. In both cases, there is only a single relevant energy scale. If 
we denote by tpp, Spg and e^^ the interaction energies for polymer-polymer, polymer- 
solvent and solvent-solvent interactions, the relevant energy scale is 

e = 2ep,-e,,-epp (11) 

The polymer-polymer interaction was set to -1 (resulting in attractive interactions for 
nearest-neighbor contacts) and the solvent-solvent and polymer-solvent interactions to 
zero. Temperature is normalized by the energy scale, so that T =kT / 1, where k is 
Boltzmann's constant. 
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For all chain lengths studied, we made certain that the box length, L, was at least 4 times 
the maximum radius of gyration of the polymers for the temperatures and densities stud- 
ied. We performed a mix of 50% particle creation/annihilation steps and 50% reptation 
steps. 

Typically, after some preliminary runs to establish the approximate location of the critical 
point, we obtained histograms for a single long run at conditions near the estimated criti- 
cal point. This was followed by additional, shorter runs near the expected liquid coexis- 
tence density at subcritical temperatures. The number of Monte Carlo steps per run was 
typically between 20 and 100 xlO^. Execution time required for a given number of steps 
is an increasing function of chain length, with a typical value being 1 hour / lO' steps on a 
hitel Pentium 100 MHz processor. The acceptance ratio for the paiticle creation / anni- 
hilation step for runs near the critical point ranged from 65 % for r=32 on the z=26 lattice 
with mean volume fraction <(p> = 0.20, to 26 % for r=1000 on the z=6 lattice with mean 
volume fraction <(p> = 0.066. These acceptance ratios are significantly higher than those 
reported by Wilding el al. ^ for the bond fluctuation model. The difference is probably 
due to the simpler character of the cubic lattices we use and the less severe volume exclu- 
sion constraints of nearby monomers. 

Data from histograms for a given lattice coordination number, chain length and system 
size were combined using the Ferrenberg-Swendsen algorithm'' and analyzed to deter- 
mine the critical point by matching to the universal order parameter distribution (see fig- 
ure 1). Coexistence densities away from the critical point, for which the liquid and vapor 
peaks of the density distributions are clearly separated, were obtained by determining the 
value of the chemical potential at a certain temperature that results in equal areas of the 
vapor and liquid peaks. The coexistence densities were determined by the first moment 
of the corresponding peak. Near the critical point, where the liquid and vapor peaks of 
the density distributions overlap, coexistence densities were established by fitting to the 
scaling relationships 
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(12) 



and 

i|)i-(|)2=tP(B + Bix'') (13) 
where 
T -T 

T = ^ (14) 

The scaling exponents corresponding to the Ising universality class, p =0.326, Ji=0.90, 
and 6 =0.54, were employed. The constants A , B , and B| were identified as the pa- 
rameters that yielded the best fit curves. The reduced distance from the critical point over 
which coexistence data were fitted to the scaling relationships was T«0. 002-0. 005 for the 
longest chain lengths and t=0.01-0.05 for the shortest chain lengths. Complete data sets 
are available on the world wide web'^. 

Statistical uncertainties for selected systems and conditions were obtained by performing 
several duplicate runs and obtaining independent estimates for the critical parameters and 
coexistence densities. 

Results and Discussion 

The first question that we needed to address is that of validation of our programs and 
analysis procedures. Some results for relatively short chains are available for the z=6 and 
z=26 lattices''^'". Our results are generally in agreement with previous literature data for 
the coexistence densities, within the stated uncertainties of the earlier results. However, 
some systematic discrepancies exist for the longer chain lengths between our data and 
those of Yan et al,' as shown in figure 2 for r=100. The coexistence liquid densities 
near the critical point are higher than the earlier results, but are probably within the sta- 
tistical uncertainty of the latter, estimated as ±0.04 in liquid volume fraction by interpo- 
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lating the error bars of figures 6 and 7 of'. Our results have statistical uncertainties 

17 

smaller than the size of the symbols used in figure 2. The data of Madden et al for the 
same system are in much better agreement with our calculations. 

Coexistence data for chain lengths r=16 to ;=1000 are shown in figures 3 and 4. Typical 
uncertainties of the calculated coexistence densities are ±1% of the reported liquid and 
vapor density. The estimated critical points as a function of chain length are shown on 
the figures by filled circles and reported in table 1. The reference point for the chemical 
potential reported in table 1 is the reversing random walk with no interactions. It should 
be emphasized that the critical point parameters are obtained for the specific simulation 
system size, L, reported in table 1. Typical statistical uncertainties for the critical pa- 
rameters reported in table 1 are ±0.001 (z=6) and ±0.01 (z=26) for the temperature, 
±0.001 for the critical volume fraction. The critical chemical potential is determined to 
within 4 significant figures for a fixed critical temperature. There is a high covariance of 
the critical chemical potential with the critical temperature, while the critical volume 
fraction is relatively insensitive to uncertainties in the critical temperature. 

The dependence of the calculated critical parameters for two specific systems on simula- 
tion system size are shown in table 2. It is clear that there is a some residual dependence 
of the critical chemical potential on system size. The dependence of the critical tem- 
perature and volume fraction on system size is below the statistical uncertainty of the cal- 
culations. According to finite-size scahng theory'" the critical temperature as a function 
of system size scales as 

r;(oo)-r;(L)«L-(«+«'^ (15) 

and the critical volume fraction scales approximately as 

((^)^(L)-(.^)^(o=)~L-0-«)'^ (16) 

Due to the presence of field-mixing in real fluids, only averages for the infinite-volume 
critical volume fraction can be obtained, as compared to direct estimate for the infinite- 
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volume critical temperature. The scaling exponents in equations 15 and 16 are 6 =0.54, 
a =0.11, and V =0.629. 

The statistical uncertainties of oiu data are too great to allow for confirmation of the fi- 
nite-size scaling functional dependence. For the purposes of analysis of the data and 
comparisons with theoretical and experimental estimates, it is sufficient to restrict our 
attention to data for a fixed system size, as reported in table 1 . 

Calculated critical temperatures are plotted as a function of chain length in figure 5. We 

have plotted IIT^ versus suggested by the Shultz-Hory relationship 

(equation 1). Linear regression of the data for r > 64 to obtain an estimate for the infi- 
nite chain length critical temperature yields ^^.(oo) =3.71+0.01 for z=6 and 
r^(oo) =20.85+0.01 for z=26. Linearity is followed very closely, with correlation coeffi- 
cients > 0.9999 for both coordination numbers. The infinite chain length critical tem- 
perature estimate is not sensitive to the "cutoff of lowest chain length included in the 
regression. For 2=6, Yan et al} report Ti.{°°) =3.45, a value 7% lower than the present 
estimate. In light of the agreement between our estimate forr^(oo) and independent es- 
timates of the temperatiu'e discussed in the following paragraph, we conclude that the 
infinite-chain length critical temperature was previously' imderestimated because of inac- 
curate critical temperatures for the longest chain lengths studied. For z=26, the present 
study is in reasonable agreement with the value T^i.^') =20 A obtained by Mackie et al. 

For the cubic lattice of coordination number z=6, Bruns'* found that two definitions of 
the temperature (Tg) are equivalent in the limit of long chains, namely (a) the second 
osmotic virial coefficient is equal to zero and (b) the mean square end-to-end distance is 
proportional to the number of bonds. The common value was obtained'* as r9= 3.713. 
Our estimate of the critical temperature in the limit of infinite chain length 
(r(,(=)=3.71+0.01) coincides with this value, thus confirming a long-standing premise 
that all three definitions of the temperature are equivalent in the limit of long chain 
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lengths. This is the first time that this important assumption of polymer solution theories 

4 

is confirmed at a level of less than 0.3% uncertainty. A previous calculation for a con- 
tinuous space model obtained agreement between the three definitions of the tempera- 
ture to within 4%. 

Scaling of critical volume fraction with chain length is shown in figure 6. For z=6, re- 
gression of the data for r>64 gives a slope of j:2=0.36±0.02, and for z=26, .V2=0.39±0.02. 
These values are in good agreement with experimental measurements'"** that give a range 
for the exponent xi between 0.38 and 0.40. These values are clearly significantly lower 
than the Hory-Huggins prediction of a:2=0.50. The only previous simulation study' to 
obtain a value for the exponent X2 comparable to the experimental value was by Wilding 
et al for the bond fluctuation model, yielding j:2=0.369. Previous studies''^ for the z=6 
and z=26 simple cubic lattices obtained significantly lower values for the exponent X2, 
primarily because they were restricted to shorter chain lengths for which the effective ex- 
ponent is lower, as can be seen in figure 6. 

Finally, a quantity of significant interest in polymer solution theories is the exponent Xi 
for scaling of the width of the coexistence curve with chain length according to equation 
3. Figure 7 shows the results of our calculations for the quantity B defined via equation 
13, as a function of chain length. The slope of the hues in figure 7 are equal to the expo- 
nent Xi. There is significant scatter in the data of figtue 7 and the exponent value depends 
on the range of data chosen for the regression. A possible explanation for the relatively 
poor data quahty is that, due to the limilations on system size that we have used, we can 
only obtain phase coexistence information up to a sizable reduced distance from the criti- 
cal point of T>0.001 (for long chains) to T>0.01 (for short chains). This may be too far for 
reMable extrapolations from the scaling relationships (equations 12-14). Despite this 
limitation, we felt it was worth analyzing the data to obtain a value for the exponent xi. 
We have chosen to use the data for chain lengths r>100 for the z=6 lattice, which yield 

xi=0.20±0.03, and the data for r>64 for the z=6 lattice, which give a:i=0.22±0.03. The 

7 19 
experimental value quoted by Dobashi et al. is ;ci=0.23. Sanchez reanalyzed the data 
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of Dobashi et al. and quoted 0.28 as the value of the exponent. Shinozaki et al. find 

20 

;i:i=0.34 from analysis of their own data. From theoretical considerations, the exponent 

X[ takes on the value 0.25 in mean-field theory, and a value of (l-P)/2=0.34 according to 

21 

de Gennes' scaling aigument. It is cleai' that additional simulation work will he re- 
quired to obtain an accurate value of the exponent X\. Such work will probably need to 
utilize larger system sizes to permit direct calculation of coexistence data in the vicinity 
of the critical point and should also investigate longer chain lengths to ensure that the 
long-chain limit has been reached. 

Conclusions 

We have used grand canonical Monte Carlo simulations combined with histogram 
reweighting techniques to obtain phase coexistence properties in polymer/solvent sys- 
tems. Simple cubic lattices of coordination number z=6 and z=26 were investigated for 
chain lengths up to f=1000. This range of chain lengths significantly exceeds the range of 
previous simulation studies for comparable systems. 

Critical temperatures were found to scale with chain length in accordance to the Shultz- 
Flory prediction. This is in agreement with previous simulation studies of continuous- 
space and lattice model homopolymers. For the z=6 lattice, the extrapolated infinite 
chain length critical temperature is 3.71±0.01, in excellent agreement with previous cal- 
culations of the temperature at which the osmotic second virial coefficient is zero and the 
mean end-to-end distance proportional to the niunber of bonds. This confirms, to an im- 
precedented level of accuracy, the standard assumption of polymer theories that all three 
definitions of the temperature are equivalent at the limit of long chains. 

Critical volume fractions were found to scale with chain length with exponent a:2=0.38, in 
excellent agreement with experimental data, but in disagreement with most polymer so- 
lution theories. Almost all previous simulations of polymer-solvent coexistence curves 
yielded lower values for this exponent because the chain lengths investigated were not in 
the scaling regime. 
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We were unable to reach definite conclusions about the scaling with chain length of the 
prefactor to the width of the coexistence curve, because of statistical uncertainties of our 
data. These uncertainties are likely the result of insufficient approach to the critical point 
due to limitations in system size. From our data, we obtain a value of the exponent xi 
near the lower range of the experimentally observed range a:i«0.23-0.34. Additional 
simulation work with larger system sizes is needed to clarify the situation with respect to 
this exponent. 

The simulation methods we have used can be used to study polymeric systems of even 
longer chain lengths at the vicinity of critical points. The reason for this is that the criti- 
cal volume fraction is lower for longer chains, thus allowing reasonable statistics for the 
grand canonical insertions and removals. The methods are also applicable to mixtures of 
different polymers or solvents. 
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Table 1 . Critical parameters as a function of chain length and lattice coordination num- 



ber. See text for estimates of the statistical uncertainties of the results. 
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0.061 


177.9 


1000 


85 


3.503 


0.072 


539.8 


85 


19.68 


0.056 


287.5 



Table 2. Critical parameters for selected system.s as a function of simulation system size. 



L 


TciL) 










z- 


=26, 


(=600 




55 


19.34 




0.068 


73.00 


65 


19.34 




0.069 


73.23 


70 


19.35 




0.068 


73.57 


95 


19.35 




0.067 


73.78 




z- 


=26, 


(=800 




75 


19.54 




0.061 


177.9 


87 


19.54 




0.062 


178.2 


110 


19.54 




0.060 


178.3 




Fig. 1 . Matching of the scaled order parameter distribution to the universal curve for the 
Ising three-dimensional universality class, indicated by the continuous line. Points are 
from our simulations for z=26: (+) r=64, L=30; (x) r=200, L=50; (O) r=600, L=95. 
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Fig. 2. Phase diagram for z=6, 100. (O) Tliis work; (■) Madden ef a/.''' 
(+) Yan et al' 



18 



0.0 0.2 0.4 0.6 

Fig. 3. Calculated phase diagrams for z=6. The estimated location of the critical points is 
given by (•) and the directly measured coexistence data by (■). Lines connect the meas- 
ured points and are extrapolated to the critical point using equations 12-14. From top to 
bottom, the curves correspond to r=1000, 800, 600, 400, 200, 100, 64, 32, 16. 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 

Fig. 4. Calculated phase diagrams for z=26. The estimated location of the critical points 
is given by (•) and the directly measured coexistence data by (■). Lines connect the 
measured points and are extrapolated to the critical point using equations 12-14. From 
top to bottom, the curves correspond to f=1000, 800, 600, 400, 200, 100, 64, 32, 16. 
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r-'/^ + 1/(2 r) 

Fig. 5. Scaling of critical temperature with chain length. Left axis, z=6: (O) This work; 
(+) Yan et al}\ Right axis, z=26: (A) This work; (■) Mackie et al} Lines are fitted to 
the critical temperatures from this work, r>64. 



